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ABSTRACT 



This thesis theoretically presents the phenomena in- 
volved in the flow of an incompressible fluid over a wedge 
with a second incompressible, lighter, and less viscous 
fluid blown through the surface of the wedge, A method is 
developed to determine the inner fluid layer thickness, the 
wall shear stress and the resulting local drag reduction. 
The results predict substantial drag reduction. 
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I. 



INTRODUCTION 



The ever increasing operating depths and speeds of modern 
submarines and the much higher resistance to the effects of 
contact or nearby underwater explosions all require a dra- 
matic increase in the speed of underwater weapons. The 
existing power plant systems seem to have reached their 
limits in power density and radical changes in body shapes 
are limited by launching facilities. The answer might lie 
in the field of boundary layer control. Possibilities are 
injection of polymers into the boundary layer to maintain 
laminar flow, suction to remove potential turbulence, more 
sophisticated weapon shapes to reduce form drag, or intro- 
ducing a gas layer between the surface of the body and the 
liquid. 

Of all the options, the two-phase boundary layer seems 
to be the most promising possibility for obtaining an order 
of magnitude increase in speed. 

In this thesis, a method is developed to calculate the 
properties of the laminar, two-phase flow produced by in- 
jecting a light, incompressible low-viscosity fluid into the 
boundary layer of the incompressible flow over a wedge. This 
is just a step on the long way to a high speed under^vater 
weapon. 
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II. BACKGROUND 



A. DISCUSSION OF RELATED PAPERS 



An extensive literature search was carried out and the 
relevant papers are discussed below. The reduction of drag 
by use of polymers is deleted from this survey because the 
drag reduction expected by this method is probably not suf- 
ficient to increase the maximum speed of the underwater 
vehicle substantially. Amazingly enough, the most recent 
paper in the field of two-fluid laminar boundary layer flow 
was published back in 1969. 

1. J. Pretsch [1] 

Pretsch starts with Prandtl's boundary layer equa- 
tions for single-phase flow with suction or blowing: 



u 



8u 

3x 



+ V 



8u 

3y 






+ V 



3^u 

^ 2 
3y 



and boundary conditions y = 0:u = 0, v = v^; yv'^ ^ : u = 0 

where u is the velocity tangential to the surface, v is the 
velocity normal to the surface, x is the arc length along 
the surface, y is the normal to the surface, U (x) is the po- 
tential velocity, v is kinematic viscositv, v (x) is the 

w 

normal velocity or blowing velocity at the surface and Re is 

Ux 

the Reynolds number. Re = — . 
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For strong suction (v^ ^ -<=°) the above equations 
simplify to a second order, ordinary differential equation 
of the form 



a^u 




+ 



du _ 



0 



y V 

where C 

V 



and boundary conditions ? = 0 ; u = 0 

r -V 00 ; u = U (x) 

The solution of this equation is the velocity profile 



u 

U(x) 



1 



e 



-C 



For strong blowing (v^ '») the assumption is made 

that the viscosity term can be neglected and the above equa- 
tions again simplify to a second order, ordinary differential 
equation but now nonlinear . 

For a number of values for similarity parameters , 
Pretsch gives analytical expressions for solutions of flat 
plate and wedge flow. 

The assumption of negligible viscosity is too re- 
strictive for the case under investigation in this thesis. 

2. E. J. Watson [2] 

This extensive report consists of three parts, i.e. 
Theory of Similar Velocity Distributions, Flow with Uniform 
Suction, and Flow with Variation of Suction Velocity. The 
report is based on the asymptotic theory of boundary layer 
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flow. Although the author is aware of the existence of 
Pretsch's work, no comparison of results is made. 

In Part I the author arrives at a number of tables 
giving velocity distribution and coefficients of the series 
expansions as functions of similarity related constants and 
dimensionless parameters. 

In Part II the similarity principle is abandoned but 
the suction velocity is uniform. The velocity is expanded 
into a series from which series for momentum and displace- 
ment thickness and skin friction are derived. In the sub- 
sequent section a number of applications are investigated 
and the results are totalized in a number of tables. 

In Part III the asymptotic theory is extended to 
cover general two dimensional flow with arbitrary distribu- 
tions of mainstream velocity and suction velocity. The 
results are again displayed in tables. The author also de- 
rives the results in Part I as a special case of Part III. 

As this report only deals with suction it is not of 
direct interest to this study. 

3. R. D. Cess and E. M. Sparrow [3] 

This paper deals with a two-phase flow in forced 
convection over a flat plate. The inner fluid is obtained 
by boiling the outer fluid at the wall so that a vapor film 
is formed. The authors arrive at the following conservation 
equations : 
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I t I 



0 



+ f f ' ' = 

g g 



0' ' + PR f 9' =0 

g 

f_ ' ' ' + f _ f ' ' =0 

Li Li Li 



where 0 = 



T - T 

00 

f 

W 00 



with boundary conditions; 

plate surface; f (0)=f '(0) = 0, 9(0) = 1 

g g 



interface : 



fL(0) 



= (1^) 
\ (PU). / 



f (n • ) 
g gi 









\(pu). / 



f " (n . ) 

g gi 



9 (n^) = 0 



free stream; f^' ' 



2 as n. 



where f^ and f^^ are the stream functions for respectively 
inner and outer fluid, T is the local temperature, is the 
temperature at the wall, is the temperature in the free 
stream, PR is the Prandtl number, p is the mass density and 
y is the dynamic viscosity. f^ is a fun ctio n of the dimen- 



sionless parameter n defined by n = ^ 



; U 






n is the 
gi 



value for this parameter at the interface of the two fluids 
f^ is a fu nctio n of the dimensionless parameter defined 



by = I 



u. 



VlX 



The subscripts L and g refer to respective- 
ly outer and inner fluid. The primes denote differentiation 
with respect to the relevant dimensionless parameter. This 
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problem is solved by using an Emmons -Leigh [4] table and a 
successive approximation procedure. 

4. K. Nickel [5] 

In this paper, the laminar flow of an incompressible 
fluid without heat transfer is investigated using a method 
of bounding the solutions applied to Prandtl's equations. 
After introducing the well-known dimensionless parameters n 
and f(n)/ Prandtl's partial differential equations reduce to 
the ordinary, third order non-linear differential equation: 

f " ' + ff ' ' + B(1 - f '^) = 0 

with boundary conditions: f(0) = f^ 

f ' (0) = 0 
f ' (oo) = 1 

The solutions are reflected in a number of graphs, 
especially "abbildung 8" giving f' (0) as a function of 3 
and fQ will prove to be very useful in this thesis. 

5. R. D. Cess and E. M. Sparrow [6] 

This paper is an extension of the paper discussed in 
Section 3 above and deals with the case of a subcooled fluid. 
To extend the earlier results it was necessary to consider 
the conservation of energy within the fluid. The results are 
again reflected in a number of figures. 

6. S. M. Sparrow, V. K. Jonsson and E. R. G. Eckert [7] 
In this paper the laminar tv/o-phase flow, obtained by 

injecting gas into a liquid at the wall of a flat plate, was 
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investigated. The authors arrive at the following ordinary- 
differential equations 



for the gas; ^g' ' ' ^g^g' ' ” ^ 

for the liquid: ~ 



with boundary conditions: 
at the wall : 



f (0) = f ^ (blowing parameter) 
g gu 



fg' (0) =0 (no slip condition) 

at the interface; f (n_^) = fr (0) = 0 (no mass transfer) 

g gx L 



f ''(n •) =Rf,''(0) (stress continuous) 
g gi L 



where R = 



(PU) 

(py) 






f ' (n .) = f_ ' (0) (velocity continuous) 
g gi L 



for n. 



f^MOy) ^ 2 (velocity approaches U ) 

1j Li 



The authors solve these equations by applying a 
series expansion and subsequently a successive approximation 
procedure, similar to the one described in the paper by Cess 
and Sparrow (Section 3 above) . 

The results are given in a table and a graph. It is 
proven that evaporation at the interface only works 
advantageous . 



7. H. P. Eichenberqer and J. D. Offutt [3] 

Under assumption of 

- shear stress in liquid is negligible 

- interfacial velocity and pressure distribution are those 
imposed by the potential flow 
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- the inertia of the gas is negligible 

- the gas flow is laminar 

- body moves steadily through water at zero incidence 

- gas and water properties are constant; in particular the 
density of the gas remains constant 

the problem of a two-phase flow over different kinds of 

bodies is solved by a finite differences method. Results 

are given in a number of graphs . 

8. W. S. Bradfield, R. 0. Barkdoll and J. T. Byrne [9] 

This paper discusses experiments with a "gassing" 

model to investigate possible drag reductions. The vapor 

layers were obtained by film boiling, sublimation of a solid 

surface and chemical reactions. The experiments indicated 

that for film boiling a region exists for which a stable 

vapor layer may be expected. 

Assuming 

- steady, two dimensional laminar flow 

- pressure everywhere constant 

- heated surface temperature constant 

- temperature of liquid vapor interface is equal to satu- 
rated liquid temperature 

- thermodynamic and transport properties of liquid and 
vapor are constant 

- thermal radiation effect on vapor and liquid layers is 
neglected 

- buoyancy effect on forced convection boundary layers is 
neglected 

- velocity and temperature profiles are linear in vapor and 
parabolic in liquid 
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and using the Pohlhausen integral method the results of the 
experiments are predicted. In some of the experiments dra- 
matic instabilities were found, and the authors had to con- 
clude that producing a stable vapor layer between a surface 
and a liquid drastically reduces friction drag but unstable, 
two-phase flow (like nucleate boiling and uncontrolled sub- 
limation) may actually cause the drag to increase. 

9. T . J . Lang [10] 

In this patent Lang "designed" a torpedo capable of 
using the gas film drag reduction phenomena. The lift is 
provided by the body shape and the exhaust gases are used to 
produce the gas layer. At the rear end of the torpedo the 
gas layer is removed by suction. 

10. E. R. van Driest [11] 

The author indicates that the quest for higher 
speeds of underwater craft can be solved by using the knowl- 
edge acquired in high speed aerodynamics . 

In a number of figures the author gives local fric- 
tion coefficients as a function of Reynolds number, speed, 
wall temperature and ocean depth. Further effects of hydro- 
dynamic heating, cavitation and insertion of additives are 
discussed . 

The author concludes that although the propulsion 
systems have to be developed, it is, in principle, possible 
to obtain high velocities for underwater vehicles. 
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B. QUO VADIS 

Considering the above discussion it was felt that more 
knowledge had to be obtained of the phenomena taking place 
in the two-phase boundary layer. A first step would be to 
extend the work on flat plates to two-dimensional bodies 
with pressure-gradients. It is logical to follow the paper 
of Sparrow, Jonsson and Eckert [7] , discussed above in Sec- 
tion A-6. 
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III. THEORY 



A. GOVERNING EQUATIONS 

To predict the skin friction and drag characterics from 
first principles the starting position is based on the basic 
conservation laws governing the velocity distribution, here 
in Cartesian coordinates, 



Continuity ; 
Navier-Stokes : 



3o , 3 ( ou) 3 ( ov) 3 ( pw) . 

at 3x 3y 3z ^ 

Du _ „ 3P ^ 3 r fo3u 2 , . ->ll 

^Dt ^ " 3x 3x[^^[^3x " 3^^^ 

, 3_ r f 3u , 3 V T , 3 f 3w 3u l 

5y]_^i3y 3xjJ 3 z|^^[Tx Tzj 

Dv _ „ 3P ^ 3 r f-3v 2 -l1 

^Dt ^ " 3y ^ 3yr[^3y ” 3^^ 



3 




3v 3w' 




, a 


* ^ 

f 3u , 3vl 


3z 


Vi 


^3z 3y, 


- 


3x 


^[3y 3xJ 



^Dw _ 3? , 3 

*^Dt ^ 3z ^ 3z 



I ~ 3w 2 , . 

u 2 ^ - rr-div W 

\ jZ J 



_3_ 




•• 

' 3w 3u) 


^ 3 


( 3v ^ 3wl 


3x 


y 


/3x 3zj 







Sqn of State : 



P = p Rg T 
y = u (P,T) 



do *^^T dW 

5^rgY_Equation: |2 = ^ 

where w is the velocity with components u, v, w in respec- 
tively the X, y and z direction, ^ is the operator 

3 3 3 D 

+ v-^ + w-^ , X, Y and Z are the body forces in 



3 Z 
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respectively x, y and z directions, P is the local pressure. 



Rg is the gas constant, Q is the heat into the system, E,p is 
the internal energy and W is the work done by the system. 

Assuming steady state, absence of body forces, viscosity 

3 P 

lj(T,P) is constant, two-dimensional flow, is negligible, 

dy 

laminar flow, and both fluids incompressible (the energy 
equation decouples) , these equations simplify to the Prandtl 
boundary layer equations: 

Continuity; 



= 0 

3x 



3y 

Navier-Stokes; u-|^ + v-|^ “ “ ^ T§ 

ox 3y p 3x „ 2 

dy 

Using the subscript g for properties of the inner fluid and 
the subscript L for the outer fluid, it is possible to write 
two sets of equations describing the two-fluid laminar 
boundary layer. 



Inner fluid: 



= 0 

3x 3y 



U |!i + = . 1 |P ^ ^ £u 

3x 3y 0 3x '^3y^ 



dy 



p 3x 



Outer fluid: |^ + |^ = 0 

3x 3y 



3x dy 



1 ^ H. 



( 1 ) 



( 2 ) 



(3) 



(4) 



To describe the problem fully, however, it is required to 
state the boundary conditions and the conditions at the 
interface between the two fluids. At the wall, in the inner 
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fluid, the velocity parallel to the wall, u, is zero (no 
slip) and the velocity perpendicular to the wall is equal to 
the imposed blowing velocity. At the interface the tangen- 
tial velocity and shear stress have to be continuous. Fur- 
thermore, the mass transfer across the interface is assumed 
to be zero. In the outer fluid u has to approach U(x) as y 
approaches infinity. A formal statement of these conditions 
is : 

at y = 0 u = 0 

V = V 

w 

at interface u = u^ 

g L 





(TO 




3Ut 1 


u 

g 


g 

3y . 


. = ’^L 
1 


L 

3y J 



m = 0 or V (interface) = v^ (interface) = 0 

g Li 



in outer fluid u U(x) as y 

For the given assumptions, equations 1-4 and the boundary 
conditions describe the problem completely. This system pre- 
sents substantial difficulties in attempting a general solu- 
tion. Similarity solutions, however, will solve the deadlock. 



B. TRANSFORMED EQUATIONS 

Similarity-type solutions are well known and have the 
characteristic that the velocity profiles have similar shapes 
at all x-positions along the surface. Mathematically this 
gives the enormous simplification that the governing partial 
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differential equations are converted into ordinary differen- 
tial equations. Before defining the dimensionless para- 
meters it should be noted that the continuity equations are 
satisfied by the stream-function ijj defined by 



u 



3x 



and V 



3x 



Following Schlichting [12] , the following dimensionless 
parameters can be introduced for the inner fluid: 



U L 



Reynolds number Re^ = — 



(where and L are some 



reference quantities) 



y/Re 



n = = — 7 —^ C - r‘ (where g(x) is a dimensionless scale 

g Lg(x) ^g L 

factor for y) 



’4j(x,y) /Re , 

f(n r — rn — \ ’Hx,y) = f(n ^ — • — 2_Ll_L 

'g'^g' L U(x) g(x) g 9 

/ q 



To simplify further manipulation, it is wise to state at 

this point that to ensure similarity f(ng,?g) can not depend 

on 5 = that f = f(n ). 

y y y 

From the above, the velocity components may be derived. 



u = u (x) f ' (n ) 
y 



and V = - 



(U(x) g(x)) + n f ' (n ) i?(x) 

/R^ /Re g g g GX 

/ q / q 



Inserting these expressions into the Prandtl boundary layer 
equation and after a considerable amount of algebra the equa- 
tion converts into 
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f + af f’’ +3(l-f^'^) = 0 

g g g 



where a = ^ g(x)) 



ft - L q (x) du (x) 

■ U dx 

00 

Reasoning in the same way for the outer fluid after 
defining 

U L 



Re^ = — 
L V, 



= (y ~ 



“ L g(x)' /^®L ' 



_ X 
’L L 



'Mx,y) 



L U(x) g(x) 



one arrives at 



f^''' + af^f^'* + 3(1 - = 0 

where a and 3 are defined the same as for the inner fluid. 

It should be noted that to satisfy the similarity re- 
quirement a and 3 should not depend on x. Still following 
Schlichting [12] , it is possible to arrive at the following 
general conditions for U(x) and g(x): 



2 - 



U(x) _ „2a-6 

“ 5 ^ ■ 



(2a - 3)^ 



2a-6 



where K is an integration constant 



g(x) = 



' U 

X ® 

(2a - 6)r 



L U(x) 



Now, without loss of generality, take a = 1 and introduce a 
new constant m = so that the expressions for and 



g(x) convert to 
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U (x) _ p,(l+m) 



U 



X 



m 



^1+m L ; 



g (x) = 



o u 

2 X «> 



V m+1 L U (x) 



Furthermore, the transformation equations for the ordinates 



are = 



/ m+1 U (x) 

2 V X 

g 



and n. 



(y-<S) 



/m+1 


U(x) 


/ 2 





With Schlichting [12], it can be noted that these re- 
quirements can be satisfied by the flow in the neighbourhood 
of the stagnation point of a wedge with included angle 3'H', 
where U(x) is given by Cx^. 

Next the boundary and interface conditions have to be 
converted by expressing them into the newly defined dimen- 
sionless parameters 



At y = 0 



or Hg = 0 



u = U(x) fg' (n) = 0 or fg' (0) = 0 



-L f (n ) , 






-V,, /Re^ 
w V g 



■g g L ^ ^ ^ 

dx 



at interface n =n- m=0orf(n-)=fT(0)=0 

g gi g gi L 



Hl = 0 



“g “ “l *g''"gi’ ' 



du 



du. 



~d^ ^g' ' <"^i^ = ^ ^ = 

for y-^-ooornT^” u-^U (x) or f,. ' 1 as Ht 

Li Li Li 



( oy) 
(py) 
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I 



I 

1 





It should be noted that 



- to ensure that independent of x; 

y y 



^ d(U(x) g(x )) 
w dx 



Inserting the expressions for U (x) 



and g(x) one obtains v,, \/ tVtv t 

^ W V U (X) 



- R depends on the properties of both fluids depending on 
temperature and pressure. 



Summarizing the problem can be formally stated by 
equations; ' ' ' 



+ f 

g 


f ’ ' + 
g 


S(1 - 


fg'"> 


= 0 


(5a) 






3(1 - 




= 0 


(5b) 


fgMO) 


= 0 










fg(0) 


= -V j 




2x 


• where v 

w 


/U(x) 


(m+1) 


v_ U(x) 


V X 






f ' ' 

a 



{T] ) = R f ’ • (0) 

y i Li 
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IV. EXPLORED WAYS OF SOLUTION 



A. SERIES SOLUTION 

Following the procedure of Sparrow, Jonsson and Eckert 
[7], a series solution of the momentum equation seems to be 
appropriate. Expanding with respect to Hg. = 0 gives 



'g'V 



f (0) + f ' (0) n + 

g g g 



2! 



+ 



f ' ' ’ (0)n„^ 

_2 2 _ + 

3 ! T . . . 



The first two terms in the expansion are available from the 
boundary conditions. The f^' ' (0) is unknown for now and will 
be determined by the solution. Higher derivatives can be 
expressed in terms of lower derivatives by taking sequen- 
tially derivatives of both sides of the momentum equation 

f ' » » = -f f - » - 3(1 - f < 2) 

and evaluating at rig = 0 . 

Some expressions for derivatives are: 



(III) 


(0) 


= -f 


(O)f' - 


(0) 


- .3 




f (IV) 


(0) 


= f2 


(0)f • • 


(0) 


+ 3f(0) 




f(V) 


(0) 


= — f 


^ (0) f ' 


' (0) 


- 2f^(0) 


+ (23-1) f* (0) 


f (VI) 


(0) 




(O)f' ' 


(0) 


+ 3f^(0) 


- (33-5)f (0)f ' 



- (63-4) Sf * ' (0) 
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It was impossible to find a recurrency in the expressions. 



On applying the interface condition = 0 makes it pos- 

sible to write 



meters. Here Sparrow, Jonsson and Eckert refer to an earlier 
paper by Cess and Sparrow [3] . This paper is discussed in 
Chapter II, Section A3. This paper gives a table relating 
f_ ' (0) and f^ ' ' (0) which corresponds to a solution of the 
equation for a flat plat. 

At this point it was realized that the data necessary to 
go any further was not available, therefore this method was 
abandoned, and other methods of solution were explored. 

B. NUMERICAL SOLUTIONS 
' 1. lODE [13] 

lODE is an interactive ordinary differential equation 
solver developed by the Computer Center of the Naval Post- 
graduate School. The program accepts a number of first order 
non-linear ordinary differential equations with associated 
initial conditions and solves the equations. 



f " (0) 

g 



-f (0) 





In this case equations 5 were converted into : 




dn ^2 



f Note; ^ 



(6a) 
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(6b) 



= y 

dn ^3 



f ' ' 



^ = £■■■ = -yiY 3 - 3 + syj" (6c 

where y^, ^2 ^3 dummy variables for the derivatives 

of f(ri)/ with initial conditions 



Yl(0) 


= f(0) 


(7a) 


Y2(0) 


= f ' (0) 


(7b) 


Y3(0) 


= f ' • (0) 


(7c) 



The equations 6 are the same for both fluids but as noted 
before the initial conditions differ. 

Sanders (14] suggested that it should be possible to 
guess the thickness of the inner fluid layer rig^ and then a 
value of f ' ' (0) could be found by trial and error so that 

g 

no mass transfer occurs across the interface (f (n •) =0). 

g gi 

The values for f , f and f' at the interface give the 

g g g 

initial conditions for the outer fluid equations 6. The 
values of f-Mny) should approach unity for large enough Ht 

1j ij ij 

while fT-' ' (Ht) approaches 0. If f_ ' (Ht) aoes not go to unity 

Ij Ij Ij Ij 

the initial guess for the layer thickness was wrong and 

another value for the inner fluid layer thickness n is 

gi 

chosen. This process is repeated until a layer thickness is 
found for which f^' (») = 1 and f^' ' (») = 0. 

The calculations carried out on the IBM 3033A? com- 
puter system showed that this procedure worked, but while the 
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calculations were fast the interactive side of the program 
has proven to be too time consuming . 

2 . Desk Top Computer 

So the necessity arose to find some way to decrease 
the time spent communicating with the computer and to auto- 
mate the trial and error process. Furthermore, it was 
necessary to find some way to interrupt the execution and 
to execute the program by steps to determine instabilities 
in the solution. 

Since a Wang 2200 desk-top computer with line printer 
was available, the problem was rewritten for this computer 
and the solution for a number of cases was carried out. It 
has to be realized that a desk-top computer is relatively 
slow so that a compromise between calculation time and ac- 
curacy of the result had to be found. It was decided that 
the primary goal was to find a method of solution and not to 
tabulate a large amount of accurate data. In the light of 
all the assumptions high accuracy would be doubtful anyway. 
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V. SOLUTIONS 



A. DEFINITION OF THE PROBLEM 

The momentum equation will be written as three first 
order differential equations as for the "lODE program": 



dn 

dn 

dys 



= y, = f’ 



= y = f” 



dn = -^1^3 - e - 6yi 
Now multiplying both sides of these equations by dn and 
changing to finite differences one arrives at 

Ayn = y^-An 



Ayo = yo-An 



Ay 3 = (-yi-Y3 - s - 3y^ ) -An 



Now by stepping n through the desired range of values, 
starting with known or guessed initial values for rig/ y^ / ^2 
and y^ it is possible to calculate y^ , ^^d y^ as a func- 

tion of n. 

It has to be noted at this point that we approximate the 
integrals by a "Lower Riemann Sum". It is known, however, 
that using a "sufficient" number of steps guarantees a satis- 
factory accuracy. This procedure is applicable for both the 
inner and the outer fluid. The initial values for the outer 
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fluid follow logically from the solution of the inner fluid 
and the interface conditions . 

Given values for the blowing parameter fg(0) and the 
properties of the fluids parameter R the procedure is as 
follows ; 

- Choose an inner fluid layer thickness. 

- Find f' (0) by an interval halving method so that f(n ■) 
(no mass transfer across the interface) equals 0 , or " 
in our program is smaller than a chosen value. 

Note; f (0) = given; f ' (0) =0 and f ' (0) is chosen as 
indicated. 

- Apply the interface conditions on the values for f, f 
and f' at the interface to find the initial values for 
the solution of the outer fluid equations. Choose an 
end value for large enough to be well "outside” the 
boundary layer. 

- Check if f^' (“) = 1 but also if ' (“) s o as it is 
possible that an unstable solution is found. In case 
the conditions are satisfied the solution is found. If 
not another guess on the inner fluid layer thickness has 
to be made. 

- The value of the inner fluid layer thickness is chosen 
in another interval halving process. 

- This process is repeated until the requirements of 

f ' (co) = 1 and f_ ' ' (“) = 0 are met. 

Li Li 

A copy of this program is enclosed as Appendix A. 

From the above, it follows that the initial intervals 

for f ' ' (0) and ri . have to be determined. The interval for 
g gi 

should be as narrow as possible and the interval for 
f ' ' (0) should be a little wider than expected to prevent in- 
stabilities and "limit-riding". These values were determined 
by trial and error and the results are discussed in a later 
chapter. 
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B. TESTING THE PROGRAM 



1 . Number of Steps of Independent Variable 

The program is used with 100 steps of the independent 
variable. A sample problem was solved with both 100 and 200 
steps. The deviations were in the fourth decimal (see Table 
5.1). The time required for the calculation however was 
doubled from approximately 45 minutes to 90 minutes. 

2 . One Fluid Flow 

It is clear that when R approaches unity, the two 
fluids are hydrodynamically identical. Computer solutions 
were obtained for a single fluid flow for 3=0 (flat plate) 
and 3=1 (stagnation flow) and further for 3 = 0.2, 0.4, 

0.6 and 0.8. The results were in good agreement with Nickel 
[5], "abbildung 8". (See Table 5.2.) 

3 . Solution for R Close to 1 

A solution with a R of 1.1 was compared to that for 
a one-fluid problem and the results were close enough to 
give further confidence in the correctness of the program to 
handle non-unity values of R. (See Table 5.3.) 

4 . Solution for Small V^edge Angles 

The solution for a 3 of 0.01 was close enough to that 
of the flat plate to give confidence about the capabilities 
of the program to correctly handle non-zero values of 8. 

(See Table 5.4.) 
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VI . RESULTS 



A. DEFINITION OF A SAMPLE PROBLEM 

To demonstrate the capabilities of the program, a special 
case was chosen for thorough investigation. For this case, 
"air at standard conditions" and water were chosen as the 
two fluids so that R = 200. For a range of blowing para- 
meters f (0), and wedge angles 6, the velocity profile was 
calculated and the layer thickness rigj_ of the inner fluid 

and the shear stress parameter f ' ' (0) at the wall were 

y 

determined . 

B. THICKNESS OF THE INNER FLUID LAYER 

Figure 6.1 gives the dimensionless thickness of the inner 
fluid layer rigj^ as a function of the blowing parameter fg(0) 
for various wedge angles 3. Figure 6.2 is an enlargement of 
these curves for small blowing-parameters. These graphs lead 
to the following observations: 

- For small values of the blowing parameter the curves are 
straight lines, ngj_ decreases very little with increasing 
wedge angle and is directly proportional to f (0). 

- For larger blowing parameters fg(0), rigj_ increases more 
rapidly with increasing fg(0). ^For large blowing para- 
meters and small angles the thickness of the inner fluid 
layer increases drastically with increased blowing 
parameter . 

- For large values of the blowing parameter, the thickness 
of the inner fluid layer decreases rapidly with increas- 
ing wedge angle 3 . 
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C. WALL SHEAR STRESS 



Figures 6.3 and S.4 give the dimensionless wall shear 
stress f ' ' (0) as a function of the blowing parameter f(0) 
and the wedge angle S . 

- For all values of wedge angle B, blowing decreases the 
wall shear stress fg' ' (0) . This reduction is spectacular 
for small values of the blowing parameter; e.g. for both 
the flat plate (6=0) and stagnation flow (6 = 1) 

f„' ' (0) is reduced to half the non-blowing value when 
fg(0) = 0.0025. 

- For all values of the blowing parameter the wall shear 
stress f ' ' (0) increases with increasing wedge angle 6; 
the smaller thickness of the inner fluid layer caused by 
the higher potential velocity. 

- For large values of the blowing parameter, the wall shear 
stress is very small and nearly equal for all wedge 
angles. The values obtained indicate that the wall shear 
stress for large blowing parameters approaches the wall 
shear stress for a wedge in an infinite volume of inner 
fluid. (See Table 6.1.) 



D. VELOCITY PROFILES 

Figure 6.5 gives velocity profiles (f'(n) vs. n) for 

various blowing parameters for a wedge angle 6 = 0.2 (36°). 

The Hr and n have been expressed in the same relative scale 

b g 

y/Ril (y-6)^/Re^ 

by noting that and resulting in 




1 

3 • 



The following observations can be made: 

- For the smaller blowing parameters the velocity profiles 
in the inner fluid are nearly straight lines. 

- The stronger the blowing the thicker the layer of the 
inner fluid causing a larger proportion of the potential 
velocity to be obtained in the inner fluid layer. For 
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small blowing parameters this proportion of the potential 
velocity increases rapidly with increasing blowing para- 
meter. Figure 6.6 gives the proportion of the potential 
velocity (or f ' (n •)) as a function of the blowing para- 
meter f (0) . ^ 

g 

- Particularly interesting is the fact that the velocity 
profiles in the inner fluid and the outer fluid have op- 
posite curvatures. This is comparable to the boundary 
layers associated with the gas-filled cavity attached to 
a body after water entry. The cavity around such a body 
has a stable interface; breakdown occurs only when too 
much gas is left at the rear separation. Although no 
rigorous proof can be given it is felt that this charac- 
teristic of the velocity profile provides stability of 
the interface in two fluid boundary layer flow. 



E. DRAG REDUCTION 



To characterize the amount of drag reduction, the local 
drag parameter fg' ' (0) obtained for the two fluid flow was 
compared to the local drag parameter obtained without 
blowing . 

To compare these local drags , a minor converson has to 
be made. The shear stress at the wall is 



T 




y=0 



In terms of the transformed variables, the expression for t 
becomes 



which can be expressed in the outer fluid parameters 



T = y 

L L g(x) ✓ L 



( py) 
(py) 



fg’Mn) 
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For the outer fluid, without blov/ing, the reference shear 
stress Tq is: 






L L g (x) / L L 



"i'>n=0 



£^( 0)=0 



Now taking the ratio; 



fg--(O) 

( 0 ) = 0 

g 



The values of this reference shear stress parameter f ' ' (0) 

1j 

are calculated and reflected in Figure 6.7. Since R is a 
known constant, for each value of the shear stress parameter 
fg' ' (0) it is possible to determine the ratio t / t ^. Table 
6.2 gives the values of t/Tq as a function of wedge angle 3 
and blowing parameter fg(0) for R = 200. 

Inspection of this table reveals that a drag reduction 
is observed for all blowing parameters and wedge angles. It 
is remarkable that for small values of the blowing parameter 
the ratio t/Tq decreases for increasing wedge angle 3 while 
for large values of the blowing parameter f(0) the ratio 
t/Tq increases with increasing wedge angle 3. For all wedge 
angles 3, however, the ratio t/Tq decreases with increasing 
blowing strength. 



F. REQUIRED POWER 

The power required to maintain blowing is 

Power = P • V 

w 

where P is the local pressure given by 
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p = p. 



dx 

Uxj 



where P 



STAG 



STAG 

is the stagnation pressure at the leading edge 



of the wedge and = -Pg m ^ < 0. If the 

inner fluid is at stagnation pressure ^ pressure dif- 

3 P 

ference of at the wall exists so that the fluid flows out 

dX 

by itself. The flow rate then can be controlled by the 
blowing hole's size. Of course, the ejected fluid has to be 



replaced requiring it to be brought up to P^ 



In practice 



STAG' 

this is of less interest as a considerable amount of exhaust 
gases of the propulsion systems at high pressures should be 
available . 
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VII. CONCLUSIONS AND RECOMMENDATIONS 



It has been shown theoretically that a dramatic drag 
reduction is possible for the infinite wedge by blowing with 
a lighter and less viscous fluid, and that the parameters 
involved can be calculated with relatively simple means. 
Furthermore the developed method gives another way of solving 
the flat plate case and the stagnation flow case. 

This work directly invites a number of follow-on investi- 
gations of phenomena associated with: 

- finite two-dimensional bodies 

- three-dimensional bodies 

- a compressible inner fluid 

- experimental verification of the theory 

- stability of the interface 
all of equal importance. 
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APPENDIX A 



"BASIC" PROGRAM FOR WANG DESK TOP COMPUTER 
TO SOLVE TWO FLUID LAMINAR FLOW CHARACTERISTICS 



1 

a 

3 

4 

5 

6 

7 

8 
S 
10 
11 
la 

13 

14 

15 

16 
17 
IS 
19 

ao 

ai 

aa 

a3 

a4 

as 

as 

a7 

as 

100 

lOS 

no 

115 



REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

REM 

DIM 

DIM 

G0( 

H6= 



THIS PROGRAM IS DESIGNED TO CALCULATE F,F',F* 

AND NGI FOR A TWO-LIQUID-LAMINAR -BOU^CARY-LAYER 
GIVEN A BLOWING -PAR AhETER F(0>,A WEDGEANGLE B AND 
A LIQUID-PROPERTIES-CONSTANT R. 

TO USE Tl-E PROGRAM ONE HAS TO GUESS ON AN INITIAL 
INTERVAL FOR THE INNER-FLUID-LAYER-THICKNESS AND 
FOR THE LOCAL SHEARSTRESS-PARAfETER F"(0). THESE 
VALUES CAN BE OBTAINED FROM GRAPHS IN THE THESIS. 
ONE HAS TO TAKE THE INTERVAL FOR THE THICKNESS AS 
SMALL AS POSSIBLE AND THE INTERVAL FOR THE SHEAR - 
STRESS A LITTLE BIT LARGER THAN EXPECTED TO PREv<ENT 
INSTABILITIES AND LIMITRIDING. 

IN SHORT INSERT VALUES FOR 

F(0) AS Gl(l), LOWER LIMIT OF F' (0> AS G6 AMD 
UPPER LIMIT OF F* (0) AS G7 IN LIHE laO 
LOWER LIMIT OF NGI AS H6 AHC UPPER LIMIT OF 
NGI AS H7 IN LINE 130 

WEEDGE ANGLE B AND FLUID PROPERTIES CONSTANT R 
IN LINE 140 

THE OUTPUT WILL CONSIST OF TWO TABLES GIVING 
FG,FG' ,FG" ,FL,FL' AND FL" AS A FUNCTION OF THE 
DIMENSI0M.es PARAMETERS NG AMD NL, PLUS A HEADING 
OF THE MOST IMPORTANT PARAMETERS AND CONSTANTS. 



GO(lOl), Gl(lOl), Ga<101), G3(101) 
L0<101), Ll(lOl), La<101), L3(101) 
1 )= 0 : Gl(l>=-.5: Ga< 1 )= 0 : G6=0: G7=a 
.95: H7=1.05 



laO G4=G&: G5=G7 
ia5 A=i: B=i: R=aoo 
130 Z=0: Z1=0 
135 H8=H6: H9=H7 
140 H=(H8+H*3)/a 
145 GS=G6: G9=G7 
150 G3<l)=<G8+G9)/a 
155 IF H=0 THEN 400 
160 IF Z < 35 THEN 180 



40 



r 



[ 




‘ 165 PRINT "AFTER 35 ITERATIONS IN GAS UNABLE TO MEET REQ'S" 
170 PRINT H,G1< 101) , GB< 101 ) , G3< 1 ) 

175 GOTO 495 
ISO H0=H/100 
185 FOR 1=1 TO 100 
190 H1=G2( I >-»1H0 

195 H2=G3( I >-»iHO 

200 H3=H0^<B-»G2( I )*G2< I >-B-A*Gl< I )-»G3( I ) > 

205 Gl< I+l )=G1< I )+Hl 

210 G2< I+l >=G2( I ) +H2 

215 G3<I+1)=G3<I)+H3 

220 G0< I+l >=G0< I ) +H0 

225 NEXT I 
230 Z=Z+1 

235 IF ABS<G1(101)> < .OOOOOOOl THEN 255 
240 IF Gl(lOl) > 0 THEN 250 

245 G8=G3(1): GOTO 150 

250 G9=G3(1): GOTO 150 

255 L0<1)=0: L1(1>=G1(101):L2<1)=G2(101>: L3< 1 ) =G3( 101 > /R 

260 M=10 

265 M0=M/100 

270 FOR 1=1 TO 100 

275 M1=L2(I)»M0 

280 M2=L3(I)*M0 

285 M3=MCr»<B-*H_2(I)*L2(I)-B-A#Ll(I)-»L3( I) ) 

290 LI ( I + l )=L1 ( D+Ml 

295 L2( I+l )=L2< I )+M2 

300 L3< I+l )=L3< I >+M3 

305 L0<I + 1)=L0< D+MO 

310 r-JEXT I 
315 Z1=Z1+1 

320 IF Zl<; 35 THEN 335 

325 PRINT "AFTER 35 ITERATI□^JS IN LIQ UNABLE TO MEET REQ'S" 
330 GOTO 415 
335 C5=L2( 101 )-l 
340 IF H=0 Tf-EN 370 

345 IF AeS<L3(101)> > .OOOl THEN 355 

350 IF ABS<C5) < . OOOOl THEN 410 

355 IF C5 > 0 THEN 365 

360 H9=H: G6=G3(1>: Z=0: GOTO 140 

365 HS=H: G7=G3(1): Z=0: GOTO 140 

370 IF AeS<L3(101)> > .0001 THEN 380 

375 IF ABS<CS) < .00001 THEN 410 

380 IF C5 > O THEN 3‘30 

385 GS=G3(1): GOTO 395 

390 G9=G3(1): GOTO 395 

395 G3< 1 )=(G8+G9)/2 

400 L3<1)=G3(1): L0(1)=0: L1<1)=GK1): L2(1)=L2(1) 

405 GOTO 260 

410 PRINT "LAMINAR FLOW OVER WEDGE WITH BLOWING" 

415 PRINT : PRINT : PRINT 
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450 PRINT "B = ",B 

455 PRINT “R = “ ,R 

430 PRINT ”FG(0) = " , G 1< 1 ) , * INI TI AL LIMITS 

435 PRINT "FGDPRIME(0)= ** , G3< 1 ) , G4,G5 
440 PRINT "NI = “ ,G0< 101 ) 

445 PRINT : PRINT : PRINT 

450 PRINT -NG" , "FG“ , "FGPRIME" , "FGDPRIME" : PRINT 

455 FOR 1=1 TO 101 STEP 10 

450 PRINT G0(I),G1(I),G5<I),G3<I) 

455 NEXT I 

470 PRINT : PRINT : PRINT 

475 PRINT "NL" , "FL", “FLPRIME", “FLDPRIME": PRINT 

480 FOR 1=1 TO 101 STEP 10 

485 PRINT L0( I ) ,L1< I ) ,L5( I ) ,L3< I ) 

490 r^XT I 

495 PRINT : PRINT : PRINT : PRINT ' •,“EM)" 
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APPENDIX B 



FIGURES 
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E’igure 6.1 Dimensionless inner fluid layer thickness as 
function of blowing parameter and wedge angle 
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small blowing 
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Figure 6.3 Dimensionless wall shear stress as a 

function of blowing parameter and wedge 
angle 
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Figure 6,4 Dimensionless wall shear stress as a 

function of blowing parameter and wedge 
angle for small blowing parameter 
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Figure 6.5 Dimensionless velocity profiles 
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Figure 6.6 Proportion of velocity obtained in the inner fluid layer 
as a function of blowing parameter 
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Figure 6.7 Dimensionless wall shear stress for the 

outer fluid without blowing as a function 
of the wedge angle 
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Table 5.2 Comparison Between Single Fluid Flow 
and Nickel (abb 8) [5] 
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Table 6.1 Comparison Between Wall Shear Stress 
for Large Bloving Parameter and Wall 
Shear Stress for Single Fluid 
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Table 6.2 Local Shear Stress Reduccion (t/T q) 
as a Function of Wedge Angle and 
Blowing Parameter for R = 200 
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